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ABSTRACT: A Thomas-Fermi- Weizsacker type theory is constructed, by means of which we are 
able to give a relatively simple proof of the stability of relativistic matter. Our procedure has the 
advantage over previous ones in that the critical value of the fine structure constant, a, is raised 
to 0.77 (recall that the critical value is known to be less than 2.72). When a = 1/137, the largest 
nuclear charge is 59 (compared to the known optimum value 87). Apart from this, our method is 
simple, for it parallels the original Lieb-Thirring proof of stability of nonrelativistic matter, and it 
adds another perspective on the subject. 
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I. INTRODUCTION 

The 'stability of nonrelativistic matter' concerns the A''-body Hamiltonian 

N 

H = Y,\Pi\ + c^V,, (1.1) 

i=l 

where Vc is the Coulomb potential of K fixed nuclei with nuclear charge Ze, with locations Rj in 
R^, and with N electrons. In units of the electron charge, e, 

V^ = -V + W + U , (1.2) 

where 

N K 

i=i j=i 

W:= Yl (1-4) 

l<i<j<N 

U:=Z^ J2 \Ri-Rj\~^- (1-5) 

l<i<j<K 

As usual p = — iV and \p\ = V—A, and the Xj are the electron coordinates. The electrons are 
assumed to have q spin states each, q = 2 being the physical value. This means that the Hilbert 
space for the A''-electron functions is the A^-fold antisymmetric tensor product of L^(R'^; C). The 
constant a = e^/^c is called the fine structure constant. 

We can easily include a magnetic field, which means replacing \pi\ by \pi + A{xi)\. The vector 
field, A, is the vector potential (in suitable units) of a magnetic field, A = curlB, and can be 
arbitrary, as far as the present work is concerned. A mass can be included as well, i.e., \pi + A{xi)\ 
can be replaced by \/\pi + ^(a;i)p + — m. The inclusion of a mass or magnetic field, while it 
changes the energy, does not affect stability. The reason for this and the requisite changes will be 
pointed out in the final section. It is for simplicity and clarity that we set m = and ^ = 0. 

'Stability of matter' means that the operator, H, is bounded below by a universal constant 
times N + K, independent of the Rj and A. In our case, because everything scales as an inverse 
length, the lower bound for H is either — oo or 0. Thus, we have to find the conditions under which 
if is a positive operator. 

Many people worked on various aspects of this problem, including J. Conlon (who gave the first 
proof [C84]), I. Daubechies, C. Fefferman, I. Herbst, T. Kato, E. Lieb, R. de la Llave, R. Weder, 
and H-T. Yau. A careful, and still current, review of the history is contained in the introduction 
to [LY88], to which we refer the reader. For present purposes it suffices to note the current state 
of affairs concerning the best available constants needed for stability, as derived in [LY88]. We can 
list these in a sequence of remarks as follows: 

1. Stability for any given values, a* and implies stability for all < a < a* and Z < Z^,. In 
fact, we can allow the nuclei to have diff'erent charges Zi, I < i < K, provided Zi < for 
all i. This follows from some simple concavity considerations and has nothing to do with the 
nature of the proof leading to a* and . 
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2. Theorem 2 of [LY88] has the strongest results, but it is limited to the case of zero magnetic 
field, ^ = 0. The result is that stability occurs if 

qa < 1/47 and Za < 2/n . (1.6) 

It is not clear to us how to incorporate a magnetic field in the proof of Theorem 2, and we 
leave this as an open problem. 

3. Theorem 1 of [LY88] has weaker results, but a simpler proof. That proof generalizes easily to 
the A 7^ CSiSG, clS pointed out in [LLS95]. The result is complicated to state in full generality, 
but a representative example is that stability holds if 

qa < 0.032 and Za < I/tt . (1.7) 

It is possible to let Za 2/7r at the expense of qa — > 0. 

4. Instability definitely occurs if Za > 2/7r, or if Z^a > 2/7r for any i. It also occurs if 

a > 128/(157r) ^ 2.72 (1.8) 

for any positive value of Z and any value of q. In other words, if a > 128/(157r) and if Z > 
then one can produce collapse with only one electron, N = 1, hy utilizing sufficiently many 
nuclei, i.e., by choosing K sufficiently large. 

5. Instability also definitely occurs if ([LY88], Theorem 4) 

a > Z^q-^I^Z-'^I^ , (1.9) 

which implies that bosonic matter (which can always be thought of as fermionic matter with 
q = N)\& always unstable. (Note: there is a typographical error in Theorem 4 of [LY88].) 



II. MAIN RESULTS 

The proof of the stability of nonrelativistic matter in [LT75] uses a series of inequalities to 
relate the ground state energy of the Hamiltonian to the Thomas-Fermi energy of the electron 
density, P{x). The chief point is the kinetic energy inequality for an A^-electron state ^, namely 

N 

(\['| X \pif I*) > const. / p^/^ . 
i=i 

The same approach will not work in the relativistic case because the corresponding inequality [D83] 
is, for dimensional reasons, 

N 

(*| X] 1*^ > const. I p'^l'^ . 
i=i ^ 

While / p^/^ can control the Coulomb attraction —Za J p(a;)/|a;|, unfortunately J p^^^ cannot do 
so. For this reason no attempt seems to have been made to imitate the proof in [LT75] of stability 
in the relativistic case. 
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However, the Coulomb singularity can be controlled by a Weizsacker type term, namely 
(Vp ) \p\ \fP)- The relativistic kinetic energy can, in turn, be bounded below by a term of this 
type plus a term of the J p'^/^ type. This and other essential inequalities will be explained more 
fully below. With the 'Coulomb tooth' now gone, TF theory with J p"^/^ can deal adequately with 
the rest of the Coulomb energy (with the aid of the exchange-correlation energy inequality [L081], 
whose remainder term also has the form J" p^/^). 

Before going into details, let us state our main results. First, we define Thomas- Fermi- 
Weizsacker (TFW) theory as follows: The class of functions ('densities') to be considered, denoted 
by C, consists of those nonnegative functions p : — > R"*" such that ^/p and ^/\p\P have finite 
L'^(R^) norms, i.e., 

C = |p : P(a;) > and J ^{1 + \p\) \V^{p)\^ dp < oo^ , (2.1) 

where \/p{p) ■= (27r)~^/^ Jj^a exp[— • x\^/p{x)dx denotes the Fourier transform of the function 
Next, we define the functional 

np) ■■= I \p\WP{p)?dp = [Vp , \p\ Vp) . (2.2) 
The TFW functional, with given positive constants (3 and 7, is then 

£{p) := f3T{p) + ^ [ p^/^{x)dx-a [ V{x)p{x)dx + aD{p ,p) + aU (2.3) 

with 

D{p ,p) := (1/2) / / p{x)p{y)\x-y\-Hxdy . 

The quantity of principal interest is the energy 

E^^^ :=M{£{p) : peC} . (2.4) 

This quantity depends on the parameters a, (3 and 7 and on the nuclear coordinates, Rj. If, 
however, we try to minimize E over all choices of the nuclear coordinates then the result is either 
or —CO, as can be easily seen from the fact that all the terms in £ scale, under dilation, as an 
inverse length. 

THEOREM 1. (Stability of TFW theory). The TFW energy, E^^^ , in (2.4) is nonnegative 
if 

(3 > -Za, and 7 > 4.8158 Z^/^a (2.5) 

On the other hand, if (3 < {Tr/2)Za then E = —00 for every choice of the nuclear coordinates. 

For the next theorem we have to define the density corresponding to an A'^-body wave function. 
If ^ is an antisymmetric function of N space-spin coordinates, normalized in the usual way, we 
define 

p.5,{x):=N / \-^{x,ai;x2,(T2;...;xN,(JN)\'^dx2---dxN- (2.6) 

l<ai,...,(7N<q 
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THEOREM 2. (TFW theory bounds quantum mechanics). Let ^' be any normalized 
antisymmetric function, with p\}, defined in (2.6). Choose 



TT 4 

B = —Za and 7 = - 
^2 '3 



1.63g-^/^ (1 - ^Za^ - 1.68 a . 
Assume that ^ is positive. Then, with this definition of the TFW functional (2.3), 



(2.7) 



(2.8) 



A corollary of these two theorems is that our Hamiltonian, H, in (1.1) is stable if 

{^)Z + 2.2159 (^I^Z'^I^ + 1.0307 (^'^ < 1/a . (2.9) 

(Cf. (1.9)) In particular, with q = 2 for electrons, relativistic matter is stable if a < 0.77 and if 
Z is not too large. When a = 1/137 the allowed Z is 59, which compares favorably with the best 
possible value 87 137(2/7r). 

We leave it as a challenge to improve our method so as to achieve the value 137(2/7r) (with 
a magnetic field present). As noted above, this value has been achieved in [LY88], but without a 
magnetic field. The most noteworthy point is the large value of the critical fine structure constant: 
acriticai > 0.77 when q = 2. 

The bound in (2.9) is, in some respects, similar to Theorem 1 in [LY88], but it is far simpler, 
clearer and gives the correct g-dependence of a (note that (1.9) gives a similar bound in the other 
direction). The chief methodological difference is that Theorem 6 is used in [LY88], which bounds 
the Coulomb potential below by a one-body potential. Here, we use the exchange-correlation 
inequality (3.9) instead. We repeat that the results above also hold with a magnetic field. 

It is to be emphasized that our stability result is really contained in Theorem 2. Theorem 
1 only gives a condition for which £{p) > 0. A better estimate on the TFW functional will, via 
Theorem 2, yield a better stability bound. 

III. SOME ESSENTIAL INEQUALITIES 

There are five known inequalities about Coulomb systems that will be needed in our proof of 
our main theorems. We begin by recalling them. 

KINETIC ENERGY LOCALIZATION, [LY88] pp. 186 and 188. 

Denote by Tj the Voronoi cell in that contains Rj, i.e., the set 

Tj := (x G R^ : \x - Rj\ < \x - Rk\ for all k} , (3.1) 

and let Dj be half the distance of the j-th nucleus to its nearest neighbor. These Fj are disjoint, 
except for their boundaries and, being the intersection of half-spaces they are convex sets. The 
ball centered at Rj with radius Dj is denoted by Bj . Obviously, Bj c Tj . 
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For any function / G L^(R^) there is the inequality 

The function Y is given by 

O 1 _|_ ^7.2 1 — 4r 

7r(l + r) 7rr(l + r2) 7rr(l + H) 7r(l + r^) 

Numerically it is found that [LY88] (2.27) 

47r / Y{rfr'^dr < 7.6245 . (3.4) 

JQ 

RELATIVISTIC KINETIC ENERGY BOUND FOR FERMIONS, [D83]. 
Let ^ and be as in (2.6). Then 

N 

{^\ V \pi\ \^) > 1.63g-i/3 / p^j3(^x)dx . (3.5) 
i=i -^R^ 



A generalization of this, of importance if we wish to include a mass, is 
{\/Pi +rn? -m\ |*) > -^m^C 

'R3 



N 

(*| [VpH^^ -"i] I*) > o^'^C / 5 ((P*(x)/C)^/3m-^) , (3.6) 
1^1 8 V / 



with C = 0.163g (sic) and with 

(7(t):=t(l+t2)V2^1 + 2t2)-^t3-ln[t + (l+t2)V2l _ (3,7) 

GENERAL KINETIG ENERGY BOUND, [C84], p.454, (and [H077] for the nonrelativistic case). 
The following bound follows from a judicious application of Schwarz's inequality. 

AT 

(*| \pi\ > (v^,b| Vp¥) • (3.8) 
This bound holds irrespective of the symmetry type of the wave function. 

EXCHANGE AND CORRELATION INEQUALITY, [L081]. If is a normalized AT-particle wave 
function there is a lower bound on the interparticle Coulomb repulsion in terms of its density: 

(*| \xi-Xj\-^\^)> D{p^,p^)- 1.68 f P^J^{x)dx . (3.9) 

l<i<j<N "^^^ 

(Once again, the antisymmetry of \E' plays no role in this inequality.) 
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ELECTROSTATIC INEQUALITY, [LY88], p.l96. First, we define a function, $ on with 
the aid of the Voronoi cells mentioned above. In the cell Tj, $ equals the electrostatic potential 
generated by all the nuclei except for the nucleus situated in Fj itself, i.e., for x in Tj 

K 

:= Z ^\x-Ri\-^ . (3.10) 



i=l 



If V is any bounded Borel measure on (not necessarily positive) then 

\( f \x-y\-^du{x)dv{y)- f ^{x)du{x) + U >\z''y^DjK (3.11) 

IV. PROOFS OF THEOREMS 1 AND 2 

To prove Theorem 1 we take /? = -kZ a/2 (if (3 > TrZa/2 we simply throw away the excess 
positive quantity). Using (3.2) with / replaced by ^/p, we have that 

£{P) > £i{p) + a£2{p) , (4.1) 
where, by adding and subtracting a term J with <E>(a;) as in (3.10), 

S^{p):=^j p'^/^{x)dx-a W{x)p{x)dx + a ^{x)p{x)dx (4.2) 

4 Jr3 Jr3 Jr3 

£2{P) := D{p ,P)- [ ^x)p{x)dx + U . (4.3) 

JR3 

The function W{x) is defined as follows: In the Voronoi cell Fj it is given by 

'Z\x-Rj\-'^, if\x-Rj\>Dj 



and 



W{x) := ^{x) + { 



(4.4) 



^ {ttZ/2)DJ^Y {\x - Rj\/Dj) , if |x - Rj \ < D, . 

Note that while the terms it J <I>p that appear in (4.2), (4.3) are merely 'strategic', the presence of 
the term $(x) in (4.4) is properly part of the potential energy of the electron and is not arbitrary. 
Actually, this strategy is the easy part of Fenchel's duality theorem (see [R70], p. 327). This duality 
principle was used by Firsov[F57] (see also [L81]) in connection with Thomas-Fermi theory; the 
full blown duality theory is not needed for our purposes, so we omit it. 

We can now seek lower bounds for £i{p) and <?2(P) separately. Using Holder's inequality, for 
example, one easily concludes that the absolute minimum of Si (p) is 

£i{P)>-f^[ [W{x)-^x)]Xdx 

47 JR3 

■f^ilY [ DT^Y{\x-Rj\/D^)Ux+ [ \x-Rj\-Ux (4.5) 
{ (f ) ' r Y{rYr^dr + 37r} Dj' (4.6) 



> - 



{azy 

{aZf 

'aZ 
4j 



*{7.6245QV3.}f:i,-'. 

9 = 1 



>-%^<! 7.6245 (^) +3n}yD-'. (4.7) 
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The last formula uses (3.4). The second integral in (4.5) is evaluated in (4.6) as Stt/Dj, and the 
explanation is the following: If we integrate \x — Rj\~'^ over the exterior of Bj we would obtain 
4iT/Dj as the result. However, we know that the Voronoi cell Tj lies on one side of the mid-plane 
defined by the nearest neighbor nucleus. This means that the integral over Tj \ Bj is bounded 
above by the integral 

DJ^ dz {27Trdr)[r^ + z^]-^ =37r/Dj . 
J I Jo 

The £2 term can be bounded using (3.11) with dv{x) = p{x)dx. Thus, 



72 K 



^2(P)>|-E^7'- (4-8) 

Combining (4.1), (4.7) and (4.8) we have proved Theorem 1. ■ 

Theorem 2 is proved by splitting the relativistic kinetic energy \p\ into (5\p\ and (1 — /?)|p|) 
with the choice (3 = -KZajl. The inequalities (3.5), (3.8) and (3.9) immediately give us Theorem 
2. ■ 

V. INCLUSION OF MASS AND MAGNETIC FIELDS 

INCLUSION OF MASS. We replace \p\ by ^Jp^ + - m and, in the corresponding TFW 
theory, we replace the right side of (3.5) by the right side of (3.6). It is not easy to carry out 
the rest of the program in closed form with this more complicated function, however. Moreover, 
it unfortunately gives a slightly worse constant than before, even when we set m = 0; instead of 
1.63q'~^/^ in (3.5) we now have C~^/^ 1.37^"^/^. The new energy will not be positive in the 
stability regime, as we had before. Instead, it will be a negative constant times N. This new value 
accords with stability and represents the binding energy of the electron-nuclear system. 

Another way to deal with the mass is to observe, simply, that \/ p'^ + m? — m > \p\ — the 
effect of which is to add a term —Nm to the energy estimate. This term satisfies the criterion for 
stability, but it has the defect that is huge in real-world terms, for it equals the rest energy of the 
electron. 

INCLUSION OF MACNETIC FIELD. Theorem 2, with a magnetic field included, is a con- 
sequence of the following two inequalities which replace (3.6) and (3.8): 

N 

(*| Y W{pi + A{xi)Y + m^-m] I*) > -m^C / g {{p^{x)/Cf'^m-^\ dx , (5.1) 
1^1 8 7r3 V ^ 



and 



N 



i=l 



To define ^J\p + A\'^ + m?' , note that if A G iL(R-^; R-^)' t^ien f ^ \\{p + A)f\\l is a closed 
quadratic form with Cq°(R^) being a form core [K7(S], [S79-1],[LS81]. Thus it defines a selfadoint 
operator and it is then possible to define ^J\p + A\^ + -m? via the spectral calculus. 
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The diamagnetic inequality for the heat kernel [S79-2] is the pointwise inequality 

I exp [-t{p + A)^] fix) I < exp [-tp''] . (5.3) 

Using the formula 

which holds for any real number a (and hence for any selfadjoint operator), we obtain the diamag- 
netic inequality for the 'relativistic heat kernel' 

exp [-iv'(pT^4)2Tm2l f{x) < exp l-t^/p"^ + m^l \f\{x) . (5.5) 

By using (5.5), and following the proof of (3.6) in [D83] step by step, we obtain (5.1). Likewise, 
(5.5) and the formula 

(/, ^{p + A)^+m^ f) = lim ^ { (/, /) - (/, exp \^-t^ {p + Ay + m^] /) } , (5.6) 

yields 

(/, \P + A\f)>{\f\, HI/I). (5.7) 
Formula (5.2) is now an immediate consequence of (5.7). 
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